! ' ' 


X-602-74-13 

. PREPRINT 


ilASA II U 


RADIATIVE TRANSITIONS 
INVOLVING THE (2p 2 ) 3 P e METASTABLE 
AUTODETACHING OF HV 


i; (<NASA~TM-X-7C570) RADIATIVE TRANSITIONS 
> INVOLVING THE <2p2) (3 Pe) KETASTABLE 
AOTOBETACHING OF H (-) (NASA) 18 p HC 
'"J$3.00 CSCL 20H 


N74-154 12 


Unclas 


G3/24 . 27232 


V. L. JACOBS 
A. K. BHATIA 
A. TEMKIN 


^ . s<: - S& 1 


M zi 


JANUARY 1974 v f. 




GSFC 


GREENBELT, MARYLAND 


) 


RADIATIVE TRANSITIONS INVOLVING THE (2p^) 3 p e 
METASTABLE AUTODETACHING STATE OF H 


V. L. JACOBS*, A. K. BHATIA, AND A. TEMKIN 
THEORETICAL STUDIES GROUP 
GODDARD SPACE FLIGHT CENTER 
NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
GREENBELT, MARYLAND 20771 


♦NATIONAL ACADEMY OF SCIENCES 
NATIONAL RESEARCH COUNCIL 
RESIDENT RESEARCH ASSOCIATE 



ABSTRACT 


The absorption coefficient for the free-bound transition H (Is) + 
e~ h” (2p 2 , 3 P e ) is calculated (together with the differential 

emission rate for the inverse process) using Is - 2s - 2p close coupling 
continuum wave functions and a Hylleraas bound state wave function. A 
maximum in the absorption and emission spectra is found to occur at a 
photon wavelength of 1219.5 A°, which is 2 A° closer to the Lyman a line 
than predicted by the calculations of Drake and is in closer agreement 
with the stellar absorption feature identified by Heap and Stecher.\rhe 
free-bound absorption process appears to be a significant source of 
continuous ultraviolet opacity. 
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I. INTRODUCTION 

The continuous photoabsorption processes 


H~ (l^S) + 

-*• H (Is) + e~ 

(1) 

H (Is) + e“ + fuo 

■+■ H (Is) + e 

(2) 


play an important role in determining the visible and infrared opacities 
of stellar atmospheres. Elaborate variational calculations (Doughty, 
Fraser, and McEachran 1966; Doughty and Fraser 1966) indicate that the 
absorption coefficients are negligible at shorter wavelengths, but the 
contributions from the autodetaching resonances of H“ which can be 
excited by ultraviolet photons were not considered. For photodetach- 
ment from the ground state, the autodetaching resonances are not 
expected to significantly exihance the absorption coefficient because 
they are the result of double-electron excitations. In the free^free 
transitions of H~ , however, the contributions from the final state 
resonances could be more significant and should be evaluated. 

2 o e 

Within the LS-coupling approximation, the doubly excited (2p ) J P 
state of H” is bound by 0.0095 ev below the H (2p) + e” continuum 
threshold (Aashamar 1969, Drake 1970, Bhatia 1970) . Because autodetach- 
ment into the H (Is) + e“ continuum would violate conservation of either 
total orbital angular momentum or parity, the (2p 2 ) 3 P e state may be 
described as a metastable autodetaching state. The close proximity of 
the (2p 2 ) 3 P e state to the excitation threshold has led Heap and Stecher 
(1974) to attribute the broad stellar absorption feature which they have 
detected on the long- wave length wing of the Lyman « line to the continuous 
(free-bound) absorption process. 
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H (Is) + e~ + H~ (2p 2 , 3 P e ) (3) 

Realistic variational calculations have been performed for the free- 
bound absorption coefficient (Drake 1974) as well as for the differential 
emission rate corresponding to the inverse (radiative decay) process 
(Drake 1973). However, the predicted spectra rise rapidly with increas- 
ing wavelength to a maximum value at 1221.5 A°, which is about 3 A° 
further from the Lyman a line than the observed absorption feature. In 
this paper, we report an independent calculation for the absorption and 
emission processes which indicates that the maximum occurs at 1219.5 A°, 
in closer agreement with the position of the observed spectral profile. 


II. THEORY 

The differential spontaneous emission rate for the (2p 2 ) 3 P e 
(Is kp) 3 P° transition, which is the only possible electric-dipole decay 
mode for the (2p 2 ) 3 P e state of H”, is given in atomic units by 

A ( U ) = l-S? 1 -§i [(Is kp,3p°|l ?! + ? 2 ||2p2,3p e )| 2 (4) 

where w is the photon energy and a, is the fine structure constant. The 
(Is kp) 3 P° scattering continuum wave function is normalized per unit 
electron energy interval . The power radiated per unit volume per unit 
solid angle and angular frequency intervals is given by the emission 


coefficient : 
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€ 


(«) = W A (w) n h~ 


(5) 


where N H - denotes the density of H“ ions in the (2p 2 ) 3 P e state. 

The effective absorption coefficient is defined in terms of the 
difference between absorption and stimulated emission and is given by 


k (w) 


4tt 3 c 2 f tfw 
i np~ [ e kt _1 . 


e (u>) 


( 6 ) 


under the conditions of local thermodynamic equilibrium (Griem 1964) . 
With use of the Saha equation 


Nh Ne~ _ 2 
n H“ 9H- 


m kt] 3 / 2 AE 
2*fi 2 } & KT , 


(7) 


where AE is the H” (2p 2 , 3 P e ) energy measured from the hydrogen ground 
state, the expression for the absorption coefficient per unit hydrogen 
atom density per unit electron density (in units of a Q 5 ) may be written 
in the form 


a (w) 


k (w) 
N H N e “ 



KT 


3/2 




A (w) 


( 8 ) 


denotes the ionization energy of ground state atomic hydrogen, and the 
electron and photon energies are related by the energy conservation 
requirement. 


E e ~ + = AE = 0.37465 a.u. 


C9) 


i 



5 


equation (8) can also be derived from the expression for the free- free 
absorption coefficient corresponding to the transition to the final state 
continuum channel H(2p) + e~ , 3 P e . Except for our inclusion of 
stimulated emission, equation (8) is equivalent to the expression for the 
free-bound absorption coefficient derived independently by Drake (1974). 

III. RESULTS AND DISCUSSION 

In our evaluation of the differential spontaneous emission rate A (w) , 
the (2p 2 ) 3 P e state of H~ was represented by an 84 - term Hylleraas bound 
state wave function (Bhatia 1970) , and the (Is kp) 3 P° scattering states 
were approximated by ls-2s-2p close coupling continuum wave functions 
(Burke and Schey 1962). The length and the velocity forms of the dipole 
matrix element in equation (4) were evaluated over the allowed electron 
energy range using the procedures developed for photoionization calcula- 
tions (Burke and McVicar 1965, Jacobs 1971) . In contrast, Drake (1973) 
has evaluated A (w) using a 50 - term Hylleraas function together with a 
discrete variational representation of and a free-electron approximation 
for the continuum. 

The position and shape of the low electron energy spectral feature 
are found to be quite insensitive to the continuum approximation but are 
found to be surprisingly sensitive to the number N of Hylleraas terms 
included in the bound state wave function. The convergence of the dif- 
ferential emission rates obtained with three different expansions is 
nevertheless satisfactory when as many as 84 terms are included (figure 1) . 
The length and velocity results presented over a more extended energy 
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range in Table 1 are in excellent agreement at low energies. For higher 
energies there is good agreement with the calculation of Drake (1973) , 
but significant discrepancies are found near the maximum which are 
described below in connection with absorption. 

The insensitivity to the continuum approximation is attributable to 
the fact that for very low electron energies the dominant contribution 
to the transition integral comes from large values of the ejected 
electron separation r where the continuum wave function may be repre- 
sented by its asymptotic form. In addition, since the triplet P-wave 
phase shift tends to zero at threshold, it was expected and found that the 
Born approximation, i.e. replacement of the ejected electron continuum 
function by j^ (kr) , gave essentially the same results as the more 
elaborate continuum functions. The same insensitivity to the continuum 
function was found by Drake (1973) . The low electron energy variation of 
A (to) is analogous to but much more rapid than the behavior of the H~ 
photodetachment absorption coefficient (Doughty, Fraser, and McEachran 
1966) , which also rises from zero at threshold to a maximum value. 

In figure 2 our 84 - term results for the free-bound absorption 
coefficient per unit hydrogen atom density per unit electron density 
are compared with the absorption coefficient calculated by Drake (1974) 
for a temperature of 5000 K°. We find a maximum at 1219.5 A° which is 
slightly larger and 2 A° closer to the Lyman a line than predicted by 
the calculation of Drake (1974) . The variation of a (to) with temper- 
ature which is illustrated in Table 2 was an important consideration in 



the identification of the (2p 2 ) 3 p e state in the stellar absorption 
spectrum (Heap and Stecher 1974) . 

Drake (1974) has pointed out that for Ne - / N H > 0.1 the free- 
bound absorption coefficient per unit length is greater than that for 
Lyman a resonance broadening (Sando, Doyle, and Dalgarno 1969) in the 
wavelength interval 1220 A° - 1300 A°. We wish to emphasize the 
dominance of the free-bound absorption process over the nonresonant 
free- free absorption process (Doughty and Fraser 1966) in the far 
ultraviolet, as illustrated by the comparison of the respective absorp- 
tion coefficients given in figure 3. 

A second maximum associated with the transition H (2p*, 5 P ) ■+ 

H~ (2s 2p, 3 P°) + hu> is found to occur at a photon wavelength of 29, 

138 A°, corresponding to radiative decay via the (2s 2p) 3 P° truly auto- 
detaching state of H”. This far infrared emission feature is shown in 
figure 4. The peak value of 4.42 x 10 3 sec” 1 a.u.” 1 is about 1/10 the 
maximum value in the far ultraviolet. Using a Hylleraas representation 
for the (2s 2p) 3 P° resonance, a procedure which does not predict the 
resonant line shape, Drake (1973) has estimated that the inclusion of 
this transition increases the integrated decay rate by only about 1%. 

This autodetaching feature would in any event be suppressed in absorption 
because of the dominance of the nonresonant free-free transitions in 
the far infrared, but it might be observable in emission under non - LTE 
excitation conditions. 

Additional sources of Lyman a resonance broadening absorption are 
expected from transitions to autodetaching resonances of H-, several of 



which are in close proximity to the n = 2 excitation threshold (Temkin 
and Walker 1965) , Such contributions to the H"" free— free absorption 
coefficient could be as significant as that found for the (2p 2 ) 3 P e 
metastable autodetaching state, and calculations for all important 
resonances are in progress . 

We are grateful to T. P. Stecher and S. Heap for bringing to our 
attention their astrophysical observation and for their encouragement 
in this investigation. We also thank Professor H, Griem for a helpful 
discussion. The use of the close coupling program written by Professor 
P. G. Burke is acknowledged. Finally, we would like to thank Mr. E. 
Sullivan for his programming assistance. 
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TABLE 1 

The differential spontaneous emission rate A (w) (in units of lO 10 
sec -1 a.u.- 1 ) for the transition H” (2p 2 , 3 P e ) -> H (Is) + e“ + tfio. The 
length and velocity results are denoted by L and V, respectively. The 


numbers in parenthesis 

are the powers 

of ten 

by which the 

entries 

to be multiplied. 





Electron energy (a.u.) 

a l 

(w) 


(o>) 

0.001 

6.576 


6.594 


0.002 

5.370 


5.395 


0.003 

4.158 


4.185 


0.004 

3.450 


3.481 


0.005 

2.981 


3.015 


0.010 

1.582 


1.617 


0.015 

9.533 

(-1) 

1.005 


0.020 

6.400 

(-1) 

6.610 

(-1) 

0.025 

4.332 

(-1) 

4.806 

(-1) 

0.030 

3.122 

(-1) 

3.740 

(-1) 

0.040 

1 . 976 

(-1) 

2.088 

(->1) 

0.050 

1.212 

(-1) 

1.279 

(-1) 

0.100 

1.727 

(-2) 

2.035 

(-2) 

0.150 

4.153 

(-3) 

5.453 

(~3) 

0.200 

1.311 

(-3) 

1.876 

(-3) 
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TABLE 2 


The free-bound absorption coefficient a (to) (in units of cm 5 ) . The 
numbers in parenthesis are the powers of ten by which the entries are to 
be multiplied. 


A (A°) 

T = 5000 K° 

T = 9000 K° 

T = 1, 

5000 K° 

1216.9 

6.49 

(-38) 

2.70 

(-38) 

1.26 

(-38) 

1217.7 

1.17 

(-37) 

4.92 

(-38) 

2.30 

(-38) 

1218.6 

1.41 

(-37) 

5.97 

(-38) 

2.80 

(-38) 

1219.4 

1.46 

(-37) 

6.24 

(-38) 

2.94 

(-38) 

1220.2 

1.42 

(-37) 

6.09 

(-38) 

2.88 

(-38) 

1221.0 

1.33 

(-37) 

5.75 

(-38) 

2.73 

(-38) 

1221.8 

1.23 

(-37) 

5.35 

(-38) 

2.55 

(-38) 

1222.7 

1.13 

(-37) 

4.95 

(-38) 

2.36 

(-38) 

1224.3 

9.60 

(-38) 

4.26 

(-38) 

2.05 

(-38) 

1225.9 

8.28 

(-38) 

3.73 

(-38) 

1.80 

(-38) 

1229.2 

6.50 

(-38) 

3.01 

(-38) 

1.48 

(-38) 

1232.6 

5.32 

(-38) 

2.53 

(-38) 

1.26 

(-38) 

1249.5 

2.13 

(-38) 

1.17 

(-38) 

6.26 

(-39) 

1266.8 

9.62 

(-39) 

6.07 

(-39) 

3.48 

(-39) 

1284.7 

4.83 

(-39) 

3.51 

(-39) 

2.15 

(-39) 

1303.1 

2.57 

(-39) 

2.14 

(-39) 

1.41 

(-39) 

1322.0 

1.41 

(-39) 

1.35 

(-39) 

9.61 

(-40) 

1361.5 

4.56 

(-40) 

5.80 

(-40) 

4.72 

(-40) 



FIGURE CAPTIONS 


Figure 1 

The differential spontaneous emission rate A (m) for the transition 
H~ (2p 2 , 3 P e ) H (Is) + e” + fttu .plotted as a function of the ejected 
electron energy and the emitted photon wavelength. N is the number 
of Hylleraas terms in the bound state wave function. 

Figure 2 

The free-bound absorption coefficient a (m) per unit hydrogen atom 
density per unit electron density for a temperature of 5000 K° , The 
results of Drake (1974) are included for comparison. 

Figure 3 

The absorption coefficients for free-bound and free-free (Doughty 
and Fraser 1966) transitions of H“ for a temperature of 5000 K°. 

Figure 4 

The resonant emission feature associated with the transition H~ 
(2p 2 , 3 P e ) -y H" (2s 2p 3 ,P°) + "Koj to the (2s 2p) 3 P° autodetaching 


resonance of H . 
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